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The ultimate climate emergency is a "runaway greenhouse" : a hot and water vapour 
rich atmosphere limits the emission of thermal radiation to space, causing runaway 
warming. Warming ceases only once the surface reaches ~1400K and emits radi- 
ation in the near-infrared, where water is not a good greenhouse gas. This would 
evaporate the entire ocean and exterminate all planetary life. Venus experienced a 
runaway greenhouse in the past, and we expect that Earth will in around 2 billion 
years as solar luminosity increases. But could we bring on such a catastrophe pre- 
maturely, by our current climate- altering activities? Here we review what is known 
about the runaway greenhouse to answer this question, describing the various lim- 
its on outgoing radiation and how climate will evolve between these. The good 
news is that almost all lines of evidence lead us to believe that is unlikely to be 
possible, even in principle, to trigger full a runaway greenhouse by addition of non- 
condensible greenhouse gases such as carbon dioxide to the atmosphere. However, 
our understanding of the dynamics, thermodynamics, radiative transfer and cloud 
physics of hot and steamy atmospheres is weak. We cannot therefore completely 
rule out the possibility that human actions might cause a transition, if not to full 
runaway, then at least to a much warmer climate state than the present one. High 
climate sensitivity might provide a warning. If we, or more likely our remote de- 
scendants, are threatened with a runaway greenhouse then geoengineering to reflect 
sunlight might be life's only hope. Injecting reflective aerosols into the stratosphere 
would be too short lived, and even sunshades in space might require excessive main- 
tenance. In the distant future, modifying Earth's orbit might provide a sustainable 
solution. The runaway greenhouse also remains relevant in planetary sciences and 
astrobiology: as extrasolar planets smaller and nearer to their stars are detected, 
some will be in a runaway greenhouse state. 

Keywords: Runaway Greenhouse, Radiation Limit, Water Vapour, Dangerous 
Climate Change, Tipping Point, Venus 

1. Introduction 

We live in a time of rapid climate change. The warming that Earth is experiencing 
now due to emissions of greenhouse gases from fossil fuel burning is unprecedented 
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in human history, and the pace of the change is very rapid compared to most 
past climate change as deduced from palaeoclimate records. It is natural to ask 
just how bad the problem could be: could humanity's collective failure to control 
our influence on climate lead to warming so cataclysmic as to render the planet 
uninhabitable? 

Wide attention has recently been drawn to such a concern by one of our most 
eminent climate scientists, James Hansen: ". . . if we burn all reserves of oil, gas, and 
coal, there's a substantial chance that we will initiate the runaway greenhouse. If we 
also burn the tar sands and tar shale, I believe the Venus syndrome [the runaway 
greenhouse] is a dead certainty" ( Hansen . 20091 p. 236). In this paper, we begin by 
reviewing the physical basis of the runaway greenhouse in order to directly address 
this issue. We then consider how geo engineering could be appropriate to address 
a threat such as a moist or runaway climate state. Lastly, we discuss the runaway 
greenhouse for other planets. 

It is a common misconception that the runaway greenhouse is a simple extension 
of the familiar water vapour feedback. As the planet warms more water evaporates. 
Water vapour is a greenhouse gas, so this enhances the warming (a positive feed- 
back). Accelerating this, by impl ication , would boil the entire o cean. Historically, 
this possibility was discussed bv lSagar] dl960h and Gold! (I 19641) . but some of the 
earliest modern-era climate models ijManabe fc Wetheraldl . 19671) showed that while 
water vapour feedback is an important positive feedback for Earth, it is not a run- 
away feedback. 

In fact, the physics of the runaway greenhouse is rather different to "ordinary" 
water vapour feedback. There exist certain limits which set the maximum amount 
of outgoing thermal (longwave) radiation which can be emitted from a moist at- 
mosphere. In the ordinary regime in which Earth resides at present, an increase 
in surface temperature causes the planet to emit more radiative energy to space, 
which cools the surface and maintains energy balance. However as a limit on the 
emission of thermal radiation is approached, the surface and lower atmosphere may 
warm, but no more radiation can escape the upper atmosphere to space. This is 
the runaway greenhouse: surface temperature will increase rapidly, finally reaching 
equilibrium again only when the surface temperature reaches around 1400 K and 
emits radiation in the near-infrared, where water vapour is not a good greenhouse 
ga s. Along the way , the entire ocean evaporates. 



Simpson! (|1927l ) was the first to i dentify a radiation li mit, finding a tropos pheric 



limit, which was seen as a paradox. Komabavashil (1967) and Ingersolll ( 19691 ) both 
analytically derived stra tospheric limits, which gav e rise to the modern concept of 
the runaway greenhouse. iRasool fc de Bergh (Il970l) applied a grey model to Venus. 
Spectra ll y reso l ved numerical calcula t ions of the full runaw ay we re performed b y 



Pollack! (l97ll) . IWatson et aB (|l984l) . lAbe fc Matsul (|l988l) and iKastind (l988). 



Over t wo decades later, these 1980s era calculations have yet to be substantially 



revised. iNakaiima et ali (jl992l ) used a grey radiative-convective model to develop 
the theoretical frame work of the various radiation limits, which is followed here. 
Pierrehumbertl (|1995l) examined the local runaway greenho use in Earth's tropics, 
highlighting the importance of dry air columns for cooling. iRennol (|1997l ) studied 
the transition towards runaway conditions (up to 330 K) with a radiative convec- 
tive model with an explicit hydrolq g ical cycle, finding mu l tiple equilibria. Pujol 



tive model with an explicit hyd rological c ycle, finding mu l tiple equilibria, rujol 
(|Puiol fc North! . 120021 IPuiol fc Fortl . 120021: IPuiol fc North) . 120031 ) and colleagues 
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elaborated a semi-grey radiative-convective model finding, for exa mple, the exis- 
tence of multiple equilibr ia during warmin g. Su givama et al. I (<2005l) examined the 
role of relative humidity. Ilshiwatari et al. ( 2002f) conducted the only study of the 
runaway greenhouse in a General Circulation Model (GC M) to date, using a sim - 
plified model wi t h the grey radiative transfer scheme from iNakaiima et al. ( 1992 ) . 



Ishiwatari et all (|2007l ) compared these GCM results with an energy balance model 



to examine climate states from low-latitude glaciation to the runaway. 



2. Physics of the Runaway Greenhouse 

In this section we review t he physics of the runa way gre enhouse in detail. Our 



theoretical framework follows lNakaiima et all (|1992h . Whilst lNakaiima et all (|l992f ) 



was very formal and mathematical, our description here is more qualitative. We 
begin by describing three different "limits" on the outgoing longwave radiation. 
Whilst we find it useful to group these together under a single umbrella term, 
we emphasize that the physical character of them is rather different. Limit 1, the 
black body upper limit is the maximum radiation that any planet of given surface 
temperature can emit; this varies with temperature. It does not lead to a runaway 
greenhouse, but is included for context. Limit 2, the moist stratosphere upper limit is 
the maximum amount of radiation which can be transferred by a moist stratosphere; 
this is invariant with temperature, so can lead to a runaway greenhouse. Limit 3, 
the moist troposphere asymptotic limit is the amount of radiation which a thick, 
pure water vapour atmosphere will emit; the radiation from any thick water-rich 
atmosphere will tend towards this asymptotic limit. It is invariant with temperature, 
so can lead to a runaway greenhouse. We then describe the changes which occur 
when increasing temperature in both runaway and moist greenhouses. 

Understanding the runaway greenhouse requires familiarity with many aspects 
of atmospheric radiative transfer and thermodynamics. Given the interdisciplinary 
interest in this topic, we summarise the necessary background material in |Appendix A| 
For those who skip |Appcndix A| most of our discussion refers to an atmosphere 
transparent to solar (shortwave) radiation and grey with respect to thermal (long- 
wave) radiation. Longwave optical depth is measured from the top of the atmo- 
sphere. Most emission to space is from where optical depth is of order unity, so 
emission depends on the temperature structure in this region. In the context of 
Earth, we refer to water vapour as a condensible greenhouse gas, carbon dioxide as 
a non-condensible greenhouse gas and radiatively inactive gases as background gas. 



(a) Radiation Limit 1: optically thin atmosphere / black body upper limit 

For an airless body (the Moon for example), or an atmosphere with no green- 
house gases or clouds, radiative emission from the surface will escape to space 
unattenuated. The flux to space at the top of the atmosphere is the black body 
emission from the surface. This flux, 

^TOA(ThinLimit) = °^surf (2-1) 

is a radiation limit: Adding any greenhouse gas (i.e. making the optical depth non- 
negligible) will cause less radiation to be emitted to space. 
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Figure 1. Radiation limits (solid lines) as a function of surface temperature, after 
iNakaiima et al\ (|l992I ). Inaccessible regions are shaded. All of the white area can, in the- 
ory be occupied depending on the amount of non-condensible greenhouse gas (e.g. carbon 
dioxide) present, but at higher temperatures outgoing flux will tend towards the tropo- 
spheric limit, as indicated by arrows. Present day mean surface temperature is 289 K with 
an outgoing thermal flux of 239 Wm -2 . Note that the change in temperature with flux 
can be seen as the climate sensitivity, so the horizontal lines certain radiation limits may 
alternatively be seen as infinite climate sensitivity and cause a runaway greenhouse. 



In Figure [T] we show this and other radiative flux limits as a function of surface 
temperature. Since the outgoing thermal flux (plotted) must balance incoming solar 
flux, and if we assume that the atmosphere is transparent to shortwave, we can 
substitute absorbed solar flux for any given planet (about 239 Wm~ 2 for Earth) 
for thermal flux and thus read surface temperature. 



(b) Radiation Limit 2: moist stratosphere upper limit (Komabayashi-Ingersoll 

limit) 

Historically, this is the limit discovered by the first workers t o explicitly describe 
the runaway greenhouse ( KomabavashiL ll967l : llngersolill969h . though as we will 
describe later, this limit is effectively impossible to achieve in practice on Earth. 
This limit is also the least intuitive. 

Consider an atmosphere with a convective troposphere and a radiative equi- 
librium stratosphere, the boundary between which is the tropopause. With the 
atmosphere in contact with an ocean, the troposphere will be moist; to calculate a 
radiation limit we assume that it is saturated with water vapour. Assuming a con- 
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Figure 2. Top of atmosphere flux with tropopause temperature. Top panel has a reverse 
scale of optical depth on y-axis, so upward on plot is increasing altitude. Bold red lines 
are solutions of the grey atmosphere model for various fluxes, showing the temperature 
profile of the stratosphere. Thin red line shows current top of atmosphere flux for com- 
parison. Blue line is optical depth of the stratosphere derived from saturated vapour 
pressure at the tropopause. Valid solutions for the tropopause are at intersections of red 
and blue lines. Bottom panel shows resulting change in top of atmosphere flux with tem- 
perature. Optical depth and grey atmosphere equations are T tr0 p = &„e s (T trop )±ff and 

ctT 4 = |-FxoA (l T + l) wnere k v is a grey absorption cross-section for water vapour, 
e s (Ttro P ) is the saturation vapour pressure at the tropopause temperature, g is accelera- 
tion due to gravity, M v / M is the ratio of the molar masses of water vapo ur and the mean 
air and T and r are an y temperature and optical depth pair on the profile l|lngersoli[l969l : 
iNakaiima et qZ.l . ll992t ). 

stant mixing ratio of stratospheric water vapour, the mass of this can be calculated 
directly from the partial pressure of water at the tropopause, which is the satu- 
ration vapour pressure for the tropopause temperature. Optical depth depends, to 
first approximation, on the mass of absorber so the optical depth of the stratosphere 
is set by the amount of water vapour it contains, which is set by the tropopause 
temperature. 

In radiative equilibrium, the upward longwave flux must equal the downward 
shortwave flux throughout. Two consequences of this are that the temperature 
increases with optical depth (this increase is rapid once optical depth r > 0.1) and 
that a higher absorbed shortwave flux will mean that the stratosphere is warmer 
throughout. 

A limit on outgoing radiation from the stratosphere arises because the radiative 



Article submitted to Royal Society 



Goldblatt & Watson 



equilibrium temperature profile needs to match the optical depth of the strato- 
sphere, which is set by the saturation vapour pressure at the tropopause. Intuitively, 
one might think first of the tropopause temperature then of the associated outgoing 
flux, but the radiation limit is easier to demonstrate by starting with the outward 
flux. For some given flux, we follow the radiative equilibrium temperature-optical 
depth profile until it intersects the temperature-optical depth curve derived from 
the saturation vapour pressure (figure [5]) . The intersection gives the tropopause . If 
the outgoing longwave flux is large, no solution can be found — the radiation limit 
is then the highest flux for which a solution can be found, around 385 Wm~ 2 here. 

One can understand this with a thought experiment in which we warm the 
tropopause. A cool tropopause means little water in the stratosphere, a small op- 
tical depth, and only a small flux is required to maintain radiative equilibrium. 
Increasing the tropopause temperature increases the amount of water vapour and 
hence the optical depth of the stratosphere, requiring a larger outgoing flux to main- 
tain equilibrium. We can think of this intuitively as a stable regime: the consequence 
of increasing temperature is an increased outgoing flux, which would maintain equi- 
librium. This regime operates whilst the stratosphere is optically thin (r < 0.1), up 
to a tropopause temperature of about 255 K. At r « 0.1, the rate of increase of ra- 
diative equilibrium temperature with respect to optical depth becomes larger than 
the rate of increase of the tropopause temperature with respect to optical depth 
(from the saturation vapour pressure) . Consequently, when tropopause temperature 
is increased, the new equilibrium has a lower outgoing longwave flux. This should 
be unstable: decreasing the outgoing flux would cause further warming, further 
decreasing the outgoing longwave radiation. Thus we have runaway temperature 
increase^. 

Thus, the simplest way to think of the Komabayashi-Ingersoll limit is to say 
that the stratosphere must be optically thin (r < 0.1), and that the radiation limit 
can be calculated assuming the tropopause is at r ~ 0.1. 

(c) Radiation Limit 3: moist troposphere asymptotic limit (Simpson-Nakajima 

limit) 

As the surface temperature rises, the amount of water vapour in the atmosphere 
increases, from water being a minor constituent (at present), to a major constituent 
and ultimately to the dominant atmospheric constituent (figure [3]). For example, 
the boiling point of water at the Earth's surface pressure of 10 5 Pa is 373 K (100°C). 
The boiling point is, by definition, the temperature at which the saturation vapour 
pressure of water equals the atmospheric pressure, so with a surface temperature of 
373 K the saturation vapour pressure of water is 10 5 Pa, the total surface pressure 
would be 2 x 10 5 Pa and the mixing ratio of water vapour 0.5 (now that surface 
pressure has increased, the boiling point would be higher than 373 K and the entire 
ocean will not boil: a liquid ocean remains until the critical point is reached at 
T c = 647 K). 

f In addition, as we move to a optically thick stratosphere, our physical assumptions break 
down. The temperature decrease with height would exceed the moist adiabatic lapse rate and 
convection would occur, causing the tropopause to move up. A very warm tropopause is unlikely 
in the real world. 
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Figure 3. Change in temperature structure (left panel) and moisture structure (right panel) 
for warming atmospheres. A background pressure of 10 J Pa is assumed, equal to that of 
Earth's atmosphere. 
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Figure 4. Temperature structure of the atmosphere with increasing surface temperature. 
A background pressure of 10 J Pa is assumed, equal to that of Earth's atmosphere. 

In general, the tropospheric temperature structure follows the moist adiabatic 
lapse rate. As the water vapour mixing ratio tends towards one, this simplifies to the 
saturation vapour pressure curve. Raising the surface temperature will evaporate 
more of the ocean — we can think of this as adding extra gas to the bottom of the 
atmosphere — but will not change the overlying atmospheric structure (figure SJ) . 

Effective emission to space takes place around t = 1. For the thick water vapour 
atmosphere we are now considering, this will be in the upper troposphere. Since the 
temperature structure is fixed, so too is the amount of radiative flux that can be 
emitted. Hence a radiation limit emerges: increasing surface temperature adds extra 
gas to the bottom of the atmosphere and increases surface pressure, but does not 
change the temperature structure of the upper troposphere from which radiation is 
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emitted to space. The radiative flux to space is decoupled from surface temperature, 
with the radiation limit 

-^TOA(TropLimit) ~ CT ^r = l (2-2) 

We refer to th is as the moist troposphere limit or Simpson-Naka jima limit, for 
Simpson ( 1927 ) who first found the limit and Nakaiima et all (1992) who described 



the relevant physics 

Calculating the value of this radiation limit requires numerical evaluation of 
a ra diative transfer code. Example values are 290 Win - 2 from a grey calcula- 
tion Nakaiima et all (|1992l ) and 310 Wm~ 2 from a spectrally resolved calculation 
(|Kastind . ll988h . 



Note that a cool atmosphere which is not pure water vapour, but contains some 
background non-condensible gas as Earth does, will "overshoot" the moist tropo- 
sphere asymptotic limit (figure [5]). This is because the presence of background gas 
decreases the water vapour mixing ratio, increasing the moist adiabatic lapse rate 
(a smaller water mixing ratio permits a bigger moist adiabatic lapse rate because 
there is more non-condensible gas to take up latent heat of condensation from the 
water) , such that the temperature around optical depth of unity is higher and more 
radiation can be emitted. Also, with a larger lapse rate but a prescribed surface tem- 
perature, the vertically averag ed atmospheric temper ature is lower, so the column 



mass of water vapour is lower (Nak aiima et all , [1992) 



Of theoretical interest, but probably not relevant to Earth, is an analogous 
dry adiabatic limit. In the preceding discussion, we have considered an atmosphere 
where the background atmosphere (10 5 Pa for Earth) is orders of magnitude smaller 
than the amount of water vapour available (~ 200 x 10 5 Pa for Earth), so the satura- 
tion water vapour mixing ratio in a warm atmosphere will always be big enough to 
make the moist adiabatic lapse rate markedly smaller than the dry adiabatic lapse 
rate. By contrast, an atmosphere with a very large reservoir of non-condensible gas 
(2 or 3 orders of magnitude more than Earth — Venus for example) would still have 
small water vapour mixing ratio in a warm atmosphere, so the moist adiabatic lapse 
rate would be no lower than the dry adiabatic lapse rate, due to the large reser- 
voir of non-condensible gas available to absorb the latent heat released. With the 
atmospheric profile tending towards the dry adiabat, the resulting radiation limit 
is hig her than the moist adi abatic limit and is almost as high as the stratospheric 
limit (jNakaiima et aLl . Il992h . 



(d) Increasing temperature 

If we imagine increasing the solar flux absorbed by the planet, the outgoing 
longwave flux must increase to compensate. How will the temperature adjust to 
provide this increased flux? The radiation limits p rovide guidance, but w e must 
interpret radiative-convective model output (from Nakaiima et al . Il992h to see 



how these limits interact (figure [5]). It turns out, perhaps unsurprisingly, that the 

f In pract i ce, th i s was the lim i t found in pre v ious numerical calculations llPollackl . Il97ll ; 
IWatson etail 1 1984 iKastinel . |198S| ; lAbe fc Matsuil . 119881 1. but was wrongly described there as 
the K omabayashi-Ingersoll limit, which should correctly refer to the stratosphere. Also. ITngersollI 
noted that a convective troposphere would not be able to deliver the required flux to meet 
the Komabayashi-Ingersoll limit, but did not sufficiently elucidate the physics. 
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Fig ure 5. Increase i n out going thermal flux as a function of surface temperature, af- 
ter iNakaiima et gl\ (1992). Black lines marked 1 to 3 show how the top of atmosphere 
flux changes with increasing surface temperature for successively higher concentrations 
of a non-condensible greenhouse gas (e.g. carbon dioxide). Line 2 corresponds to Earth's 
present amount of non-condensible greenhouse gases and lines 1 and 3 are illustrative 
of lower and higher concentrations, respectively. All lines are for an amount of back- 
ground, non-condensib le and non-radiatively active gas similar to Earth's (see figure 6 of 
INakaiima et all (|l992h for other background gas inventorie s). Radiation limits are shown 
in colour (see figure [5] for labels). Figure 9 of Kasting ( 1988j), derived from a spectrally-re- 
solved model, has similar features. 

path surface temperature takes depends strongly on the amount of non-condensible 
greenhouse gas, such as carbon dioxide, which is present. 

If the atmosphere has a negligible amount of carbon dioxide, then with a low 
surface temperature the atmosphere is optically thin and follows the black body 
limit. The temperature and outgoing flux increase together in this way whilst the 
atmosphere is not optically thick, overshooting the moist troposphere limit. The 
flux approaches, but never reaches, the stratospheric limit (or dry troposphcric 
limit). As the atmosphere becomes warmer and richer in water vapour, the lapse 
rate decreases in magnitude, meaning that it is cooler around r = 1 and less flux 
is emitted to space. The outgoing flux then asymptotes to the moist troposphere 
limit (figure [5J . Clearly, decreasing outgoing flux with increasing temperature is an 
unstable situation, so the onset of the runaway greenhouse is at the turning point 
of the flux with respect to temperature. 

Now consider an increased carbon dioxide abundance. With a fixed absorbed 
solar flux, more CO2 gives a warmer surface; this is the familiar greenhouse effect. 
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Wavelength (/Um) 

Figure 6. Absorption spectrum of water vapour 0.3 - 40 microns, shown at 220Pa and 
260K. Note the 'window' regions where the absorption coefficient is low and the general 
decline of absorption coefficient at shorter wavelength. 



If more CO2 is present whilst we increase surface temperature the flux overshoots 
the moist troposphere limit less: the runaway greenhouse is initiated with a lower 
flux, but initiation occurs with a higher surface temperature. With a suitably large 
CO2 inventory, there is no overshoot and the moist troposphere limit is approached 
monotonically (figure [5]). In no case is the runaway initiated with a flux lower than 
the moist troposphere (Simpson-Nakajima) limit. 

In the simplified system, where the atmosphere is transparent to solar radiation 
and grey with respect to thermal radiation, the runaway greenhouse occurs when 
increasing surface temperature is accompanied by constant or decreasing outgoing 
longwave flux. This implies that the surface should warm infinitely There is an 
end however, and temperature will stabilise again. To explain this, we need to re- 
lax our assumption that water vapour is a "grey" absorber: in reality it is quite 
far from being so. The absorption spectrum of water vapour (figure [6]) shows that 
there are major "window" regions between 8 and 12 fim and around 4 /im, where 
absorption is weaker than elsewhere, so radiation can escape from the surface or 
lower atmosphere directly to space. Today, the wavelength of maximum emission 
is 10 /im and much of the radiation that escapes Earth does so through the 8 to 
12 /im window. With more water vapour in the atmosphere, the 8 to 12 /im window 
will becomes optically thicker, absorbing more radiation and preventing escape to 
space. However, as temperature increases, the wavelength of maximum emission 
decreases in inverse proportion (Wien's displacement law), and radiation can es- 
cape through the shorter wavelength window regions. This allows escape from the 
runaway greenhouse. Numerical calculations shows that this occurs when the sur- 
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face is around 1400 K (Kasting, 1988; Ab e fc Matsui . 19881) when the corresponding 
wavelength of maximum emission would be 5 /im. Making the analogy to our orig- 
inal grey case, allowing emission through the 4 /im window can be thought of as 
analogous to reducing the grey absorption coefficient of the atmosphere when the 
surface temperature exceeds 1400 K. 

1400 K is well above the critical temperature of water, 647 K, at which liquid 
and vapour phases cease to be distinct, so the oceans all evaporate during the 
runaway greenhouse (the exit temperature from the runaway is basically fixed by 
the position of water vapour windows, so does not depend much on the volume of 
the oceans or whether they have all evaporated). 

Relaxing the assumption that the atmosphere is transparent to shortwave radi- 
ation causes some qualitative changes as water vapour absorbs solar radiation. As 
the water vapour inventory increases with temperature, there is initially an increase 
in absorption (decrease in planetary albedo). However, as the atmosphere becomes 
thicker, increased Rayleigh (molecular) scattering dominates and net absorption de- 
creases (planetary albedo increases) . Consequently, the runaway greenhouse would 
be initiated later, as solar absorption is typically decreasing as thermal emission 
asymptotes ( Kastinel 19881) . 



(e) The moist greenhouse 

So far, we have considered only the strict case of a runaway greenhouse. However, 
there exi s ts an intermediate atmospheric state, the so called "moist greenhouse" 
( Kastinel 1988) in which water vapour is a major constituent of the troposphere 



and the stratosphere becomes moist (figure[3]). For the troposphere, the explanation 
is simply the increase in saturation vapour pressure with temperature: with a 320 K 
surface, for example, the water vapour mixing ratio at the surface would be 10%. 
Understanding the stratosphere requires slightly more physics. We can assume that 
the stratospheric water vapour mixing ratio is determined by the mixing ratio 
at the tropopause. This is known as the cold trap; the minimum in temperature 
causes a minimum in satruation vapour pressure, which keeps water out of the 
stratosphere. However, with a warmer surface, the tropopause moves upward and 
to a lower pressure, but stays at a fixed temperature. Water vapour mixing ratio 
is the ratio of the vapour pressure to ambient pressure, so with a lower ambient 
pressure the tropopause mixing ratio becomes higher and the cold trap less effective. 
Increasing surface temperature from 290 K to 300 K, 310 K or 330 K gives increases 
in stratospheric water mixing ratio of factors of 1.5, 2.7 and 21 respectively. 

In principle, a moist greenhouse is a stable climate — just one much hotter than 
today. This is not a cause for comfort however, as the differences relative to today 
would be profound. With a warming climate, the latitudinal temperature gradient 
would weaken (likely melting any year round ice) and the atmos pheric circulation 



patte rns would be quite different (compare figures 8 and 9 of llshiwatari et at 



20021 ). Furthermore, the present stratosphere begins at around 100 hPa and the 
ozone maximum is around lOhPa. In our simple model (figure [3]) which tends to 
put the tropopause too low, by the time the surface temperature is 330 K the 
tropopause height would have risen to engulf the ozone maximum. This could de- 
stroy th e ozone layer via catalytic ozone destr uction by water derived hydrogen 



radicals ( Kasting fc Donahue . 19801 Larv . 1997 ). Even before this, much more wa- 
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ter in the stratosphere due to a weaker cold trap will rather change the stratospheric 
chemistry. A moist greenhouse would also increase climate s ensitivity. Mor e water 
in the stratosphere will increase the rate of hydrogen escape (|Huntenl . fl97^ which, 
although not relevant on human timescales, is a major process in planetary evolu- 
tion. 

Getting to a moist greenhouse simply requires warming, n ot exceeding radiation 
limits : more greenhouse gases could certainly get us there. iKasting fe Ackerman 
(|1986f ) found that somewhat over 10,000 ppmv CO2 was required. 



3. Could anthropogenic global change cause a moist or 
runaway greenhouse? 

(a) Runaway greenhouse 

Figure [3] illustrates how increasing the carbon dioxide inventory of the atmo- 
sphere affects the change in outgoing longwave radiation with temperature, using 
a grey atmosphere model. Presently, Earth's absorbed solar flux is smaller than all 
of the radiation limits described above and the surface temperature adjusts so that 
outgoing longwave radiation matches the absorbed solar flux. More carbon dioxide 
means that the surface must be warmer to provide the same outgoing flux — this is 
the familiar greenhouse effect. 

The runaway greenhouse only occurs when the outgoing longwave flux reaches 
a radiation limit. The fundamental point is that adding carbon dioxide does not 
increase the outgoing longwave flux, so cannot cause a runaway greenhouse. Whilst 
this result comes from simple models, a qualitatively similar result can be obtaine d 
from spectrally resolved radiative transfer codes (see figure 9 of Kasting . Jl988): 



even 100 bar o f CO2 does not give a radiation limit ( Kasting fc Ackermanl 1 198(1 
Kastinel[l988h . 



And yet, because of the extreme seriousness of the consequences, this theoret- 
ical answer is insufficiently satisfying. Whilst we can construct simple models to 
describe the limits of how the atmosphere behaves (and verify them with numerical 
calculations to an extent), we actually have little idea of exactly how the atmo- 
sphere will change when subject to extreme heating. Is there any missed physics 
or weak assumptions that have been made, which if corrected could mean that the 
runaway is a greater risk? We cannot answer this with the confidence which would 
make us feel comfortable. 

One major uncertainty is relative humidity. We have assumed a saturated tro- 
posphere as the end member which makes a runaway greenhouse most likely, with 
the aim that uncertainty here would be how much less likely the runaway green- 
house would be. If Earth's tropics were saturated, they would be in a local runaway 
greenhouse, whereas in reality columns of unsaturated air allow radiation to escape 



(|Pierrehumbertl I1995I ). How relative humidity will change with warming is very 
poorly understood, but will have a first order effect on temperature change. Ex- 
plicit inclusion of a hydrological cycle in models has given rise to multiple equilibria 



riRennd . fl997h . which adds an additional threat. 



Likewise, we have not discussed clouds (the physics is hard enough without 
the m). Presently, the g reenhouse effect of clouds is ab out half the al bedo effect 



(e.g. Zhang et ai . 2004). A typical argument, following IKasting (1988), is that in 
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an optically thick atmosphere, the cloud greenhouse effect would diminish as non- 
cloud optical depth would be so high and the albedo effect would dominate further, 
making the runaway greenhouse less likely. This is not a watertight argument: for 
example, with much more water vapour aloft high clouds might become thicker and 
more widespread, causing net warming. We do not know what will happen. 

There are uncertainties in the r adiative transfer too. The last full, s pectrally 
resolved, treatments of the problem ( Abe fc Matsui , 1988 ; Kastine , 19881) are over 
two decades old. Much of the uncertainty relates to the so-called continuum absorp- 
tion of water vapour (in window regions, both 8 to 12 /j,m and 4 /zm, there is more 
absorption than one would expect considering nearby water vapour lines). This is 
probably caused by far-wing line absorption by water lines, which depends on self- 
collisions of water molecules. Hence this depends on vapour pressure squared, and 
is highly relevent to our problem. Understanding and empirical constraints on this 
continuum have improved, but remain unsatisfactory. 

Are there alternate ways that the radiation limit could be reached? Could de- 
creased albedo increase the absorbed solar radiation sufficiently? Taking the lowest 
radiation limit to be 285 Wm~ 2 , one would need to reduce the planetary albedo 
from 0.3 to 0.16 to bring the absorbed solar radiation up to the radiation limit. In 
radiative forcing terms, this is 45 Wm -2 . C hanging all land surface to ocean (which 
is darker) would give a forcing of 5 Wm~ 2 ( Goldblatt fc ZahnTel 2011 ). Clouds are 
more important: removing all low cloud would give a forcing of 25 Wm -2 . Marine 
stratus depends strongly on boundary layer processes which is are not represented 
well in all models. Whilst albedo changes alone would unlikely be sufficient to cause 
a runaway greenhouse, a major reduction in cloud reflection would pose a hazard. 
In a warmer climate, more atmospheric water vapour would increase the amount 
of incoming solar radiation absorbed in the atmosphere, also decreasing planetary 
albedo. 

As a second alternative, we consider whether anything might lower the radiation 
limit. Given that the runaw ay greenhou se threshold depends o n absorption in the 
8 to 12/xm ( Kastine . 1988) and 4^m ( Pierrehumbert . l2010h windows, it seems 
plausible that a greenhouse gas which absorbed strongly here (carbon dioxide does 
not) could lower the radiation limit. This has not been investigated. Alternatively, 
more extensive high cloud cover could block radiation through the windows (clouds 
do approximate well to a grey absorber). If such clouds had larger particle size than 
present, their albedo could remain low. 

Although it seems quite unlikely, a case of more widespread and larger droplet 
size high clouds or good gaseous absorbers in the water vapour window and a 
reduction in low clouds could, in principle, increase the absorbed solar radiation and 
decrease the radiation limit. Without better understanding whether such changes 
are possible, we cannot totally rule out a runaway greenhouse. 



(b) Moist greenhouse 

The above considerations all apply to the case of a strict runaway greenhouse. 
A transition to a moist greenhouse (*Qfej) or other hot climate state is not excluded 
by theory and must be seen as a potential threat until proved otherwise. In our 
under standing , this is the physically correct interpretation of the severe hazard of 
which Hanse rJ (|2009h warns. 
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The question here is simply how much could hum an action increase the strength 
of the greenhouse effect? Kasting fe Ackerman ( 19861) found that, with carbon diox- 
ide as the only non-condensible greenhouse gas, over 10,000 ppmv would be needed 
to induce a moist greenhouse. This is likely higher than could be achieved than 
by burning all the "conventional" fossil fuel reserves — though the actual amount 
of fossil fuel available is poorly constrained, especially when one includes "exotic" 
sources such as tar sands (which are already being exploited). Greenhouse gases 
other than carbon dioxide, cloud or albedo changes could all contribute further 
warming. L ikewise , the exhibition of multiple equilibria in the relevant tempera- 
ture range (jRennd . Il997t iPuiol fc Nor 2002) complicates matters. 



(c) Could we predict this? 

If we did approach a moist or runaway greenhouse, how would we know? If 
the approach to a radiation limit is monotonic, then increased climate sensitivity 
might serve as a warning: close to a limit, any forcing (e.g. from more CO2 ) would 
cause a much enhanced temperature change. To put bounds on this, models used 
for the IPCC reports generally give a climate sensitivity to doubling CO2 of around 
3K with a likely range of 1.5 to 4.5 K. Ensemble experiments in which uncertain 
parameters are varied over reason able ranges yield a dist ribution of sensitivities 



with a long tail out to over 10 K ([Stainforth et all 120051 ). Such high sensitivity 



is inconsistent with our knowledge of palaeoclimate and the model cases which 
provide the extremes do not seem likely (due to poor representation of contemporary 
climate): thus if such a sensitivity was observed, it would indicate moving to a 



markedly different climate state. However, figurelBlffollowing lNakaiima et cd. 



and various other stud i es of the approach to the runaway greenhouse " jRe 




llshiwatari eTail , l2002t IPuiol fc Northl l2002t ISugivama etal\ . 120051) indicate that 
there are strong non-linearities and some bifurcation points (transitions between 
different steady states) and multiple equilibria on this approach. The approach to 
a runaway greenhouse is not sufficiently understood, so it is not obvious what we 
should look for as warning. 

Ideally, we would want numerical climate models to robustly resolve the transi- 
tion to a much hotter atmosphere. However, most such models have been developed 
for fairly small perturbations from the existing climate and their wider applicability 
may be limited by the obvious unavailability of data to tune the model to, and by 
simplifications made to reduce computational cost (for example, the performance 
of the radiative transfer schemes used in some contemporary general circulation 
models deterio rates at or soon after doubling carbon di oxide from present day 
concentrations ( Collins et all 120061: iGoldblatt et all 120091) ) . Considering an atmo- 
sphere that substantial fractions of which can condense introduces many additional 
challenges. A new generation of model may well be needed. 



4. Geoengineering solutions to a Runaway Greenhouse 

Suppose that the runaway greenhouse did prove to be a risk, either in the fore- 
seeable future for humanity (if our analysis above is erroneous) or for our remote 
descendants as solar output increases. Geoengineering would be necessary to keep 
Earth habitable, but how would this best be done? Three principles might guide 
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action: (1) The threat of a runaway greenhouse would come from external heating 
(probably solar), not from greenhouse gases. (2) This would be a permanent threat. 
(3) Climate sensitivity would be very high. 

From the first point, it is clear that the geoengineering strategy would be to re- 
duce the amount of sunlight absorbed. From the latter two points, it is equally clear 
that this strategy would need to be long lived and robust. Temporary reflectance- 
altering strategies such as addition of aerosols to the stratosphere might provide 
a quick fix in dire emergency, but would not answer the need for permanence. A 
more robust approach would be to deplo y sunshades at the inner Sun-Earth La- 
grange point (LI). A scheme described by Angel (2006) would use millions of small 
sunshades at the LI point, which would block 1.8% of the incident solar radia- 
tion (equivalent to about 6Wm~ 2 when averaged over the Earth). However, the 
small "flyers" would also have a limited lifetime, probably only measured in decades 
(|Angeil2006h . so constant maintainance would be required. 

If the risk of a runaway greenhouse is distant in time, we can consider solu- 
tions beyond our present engineering capability. To counter the brightening Sun, 
one could move the Earth's orbit outward (astronomical engineering). Physically, 
this could be done by altering the orbit of an asteroid to perform gravity as- 
sists around Earth and J upiter, transferring orbital energy from Jupiter to Earth 
( Korvcanskv et all . l200lh . This would be ideal as it would directly address the cause 
of the problem and would not require continual maintenance. 



5. The runaway greenhouse on other planets and exoplanets 

The theory relating to the runaway greenhouse was developed in the planetary 
sciences, concerning the evolution of Venus and planetary habitability. Reaching a 
radiation limit is seen as setting a hard limit on the inner edge of the circumstellar 
"habitable zone" (the region in which an terrestrial planet w ould likely to be hab- 
itable for life as we know it ( Abe . 1993 ; Kasting et al , 19931 )) and is thus relevant 
for the design of exoplanet characterisation missions. 

A runaway greenhouse is caused by an increase in external heating. One heat- 
ing source is solar: stellar luminosity increases with age (e.g. the sun is 30% more 
luminous now than when it entered the main sequence), so time is sufficient. Alter- 
nately, if the orbit of a planet evolves inward, a radiation limit could equivilently 
be reached. Theoretically, a decrease in albedo could mimic inc reased solar flux 



(e.g. deglaciation from total ice cover whilst solar flux is high (jlshiwatari et all . 
120071)). The second comm only discussed heating sorce is a large astroid impact 



(e.g. lAbe fc Matsui . 19881) . Impact heating is effectively delivered instantly, giving 
an immediate transition to runaway conditions; these would be transient with the 
atmosphere slowly condensing afterwards. 

For insolation driven runaways (or a planet subject to relentless bolide impact), 
the ultimate escape from a runaway greenhouse would be loss of water to space 
by disassociation due to photolysis and then escape of hydrogen to space. Given 
that the atmosphere would be dominated by water vapour, hydrogen escape would 
be energetically limited (hydrodynamic escape) and the timescale for loss of an 



Earth -size ocean of water would be a few hundred million years (|Watson et al 
19811) . 
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These considerations are not entirely hypothetical: Venus' atmospheric state 
is thought to be explained by a runaway greenhouse and subsequent hydrogen 
escape. Its present atmosphere is 92 bars, consisting 96.5% CO2 and 3.5% N 2 . 
Water vapour is present at only ~ 30ppmv ( Meadows fc Crisp . 19961) . Being nearer 
to the sun, Venus recieves nearly twice the energy from the Sun as Earth does, so 
will have exceeded the radiation limits at some point in its history. A hundred-fold 
enrichment of the deterium-hydrog en ratio on Venus rela tive to Earth is evidence 
of the subsequent hydrogen escape ( Donahue et ai . 19821) . 



The radiation limits vary with the acceleration due to gravity, g (bigger for heav- 
ier planets). We focus here on the moist troposphere asymptotic limit (Simpson- 
Nakajima limit), as it is the lowest and most commonly encountered limit. To first 
order approximation, the optical depth of a layer depends on the mass of absorbing 
gas in that layer. With this approximation, the required mass of absorber in the 
layer between the top of the atmosphere (r = 0) and the effective emission level 
(r = 1) is constant. If g is increased, the pressure at the bottom of this layer (at 
t = 1) will similarly increase (p = force/area = mass • g /area). With the vertical 
temperature structure following the saturation vapour pressure curve, higher pres- 
sure at the effective emission level corresponds to a higher emission temperature, 
hence a higher limiting ntuU A Mars-size planet would have a lower limiting flux 
and be more susceptible to a runaway greenhouse, whereas a few Earth-mass planet 
would have a higher limiting flux and be less susceptible to a runaway greenhouse^] 



6. Conclusion 

The runaway greenhouse is well defined in theory: As the surface warms, the at- 
mosphere becomes optically thick with water vapour, which limits the amount of 
thermal radiation which can be emitted to space. If the planet absorbs more solar 
energy than this limit, runaway warming ensues. The practical limit on outgoing 
thermal radiation occurs when the atmospheric structure tends towards the satura- 
tion vapour pressure curve of water as the atmospheric composition tends towards 
pure water vapour, giving a limit of around 300 Wm~ 2 . Earth presently absorbs 
around 240 Wm~ 2 of solar radiation. Increasing carbon dioxide concentration will 
make surface warmer with the same outgoing thermal flux. Following this theory, 
we are not near the threshold of a runaway greenhouse. However, the behaviour of 
hot, water vapour rich, atmospheres is poorly understood and much more study of 
these is necessary. 

In the event that our analysis is wrong, we would be left with the situation in 
which only geoengineering could save us. The only useful methods would be those 
which would reduce the amount of solar energy that the Earth absorbs. 

f Optical depth does actually depend on pressure too due to "pressure broadening" of absorp- 
tion lines. Thus, in a more complete trea tment, the limiting flux varies more slowly with g than 
indicated here. See iPie rrchumbcrt (2010) for a quantitative treatment of this. 

X A semi-analy t ical t reatment for the moist troposphere limit has recently been developed 
by IPier rchumbcrt (20l3), by substituting the vertical structure given by the saturation vapour 
pressure curve into a two-stream radiative transfer equation. Coefficients for radiative transfer 
terms are tuned to runs of a radiative transfer code. This formulation will likely be useful for 
ad-hoc calculations for applications of the problem (for example calculating a bound on the inner 
edge of the circumstcllar habitable zone for variable size extra-solar planets). 
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Lastly, we emphasise that nothing in this article should detract from the fact 
that anthropogenic greenhouse gas emissions are a major threat to global society 
and ecology. Whilst we make the technical point that these emissions are unlikely 
to cause a "runaway greenhouse" , very severe and dangerous climate change is a 
real possibility. Abrupt change to a hot "moist greenhouse" state, which though not 
technically a "runaway" would nonetheless be dire, is not excluded. The imperative 
to cut greenhouse gas emissions remains. 

Thanks to Kevin Zahnle and Ray Pierrehumbert for enlightening discussions on the run- 
away greenhouse, to Kevin Zahnle, Naomi Vaughan, Megan Smith and Tyler Robinson for 
comments on the manuscript and to Jim Kasting and an anonymous reviewer for construc- 
tive reviews. Thanks to Richard Freedman for plotting figure [5] CG was funded by NASA 
Planetary Atmospheres program (grant NNX11AC95G) and through the Virtual Plane- 
tary Laboratory Lead team of the NASA Astrobiology Institute. AJW was supported by 
a Royal Society Research Professorship. 



Appendix A. Technical Background 

In this appendix, we summarise the background physics required to understand 
the runaway greenhouse, for the benefit non-experts. Key terms are emphasised for 
eas y reference. Furthe r back ground can be found in various textbooks. For ex am- 
ple, IWallace fc Hobbd (|2006h at a n introductory level. ICurrv fc Webster] (jl999h for 



reference on thermodynamics and IPierrehumbertl (|201fJ for an accessible, interme 



diate level, coverage of radiative transfer and discussion of the runaway greenhouse. 

We are concerned with how the atmosphere absorbs and emits radiation. Two 
basic properties here are the absorptivity, A, and emissivity , e. Absorptivity is 
the fraction of the radiation at any given wavelength which is absorbed by some 
atmospheric layer and emissivity is the fraction of the maximum radiation that the 
layer could emit (the black body radiation, to be discussed shortly) which is actually 
emitted. For any given wavelength, absorptivity equals emissivity (Kirchhoff's 
law). A body which is a perfect absorber and emitter (A = e = 1) at all wavelengths 
is described as a black body. The wavelength distribution of emitted radiation from 
a black body depends solely on the temperature of the body. This distribution is 
called the Planck function. Two properties of this function are most relevant 
to the atmosphere. First, the integrated radiative emission (the radiative flux) is 
given F = aT A (the Stefan- Boltzmann law), where a = 5.67 x 10~ 8 is the Stefan- 
Boltzmann constant and T is the temperature. Thus the emitted radiative flux is a 
very strong function of temperature. Second, the peak in the Planck function (the 
wavelength of maximum emission) is inversely proportional to the temperature of 
the body (Wien's displacement law). 

The Earth's surface and the Sun correspond roughly to black bodies with tem- 
peratures of 289 K and 5900 K, which have their respective Wein peaks at 10 /im 
(thermal infrared radiation) and 0.5 //m (visible radiation). Hence terrestrial radi- 
ation is often referred to as thermal or longwave radiation, and solar radiation 
as shortwave radiation. It is common to treat these two wavelength regions sep- 
arately as their behaviour in the atmosphere is rather different. In the shortwave 
region, reflection dominates over absorption. Around 30% of incident solar radia- 
tion is reflected by clouds, the surface and by molecules (Rayleigh scattering) 
but only a small fraction is absorbed. Neglecting this absorption is a common ap- 
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proximation. In the longwave, scattering (reflection) is negligible but absorption is 
very important. This absorption is done both by clouds and by gases (greenhouse 
gases). Whilst radiative absorption has a very rich wavelength dependence, a lot of 
theoretical progress can be made with the simplifying assumption that each wave- 
length region (long- and shortwave) can be described with a single emissivity. This 
is the grey atmosphere approximation, and is used for most of the discussion in 
this paper. 

Now, consider the energy balance of Earth. Absorbed solar radiation must equal 
emitted thermal radiation at the top of the atmosphere, so we can write j(l — a) = 
aT*, where S = 1368 Wm~ 2 is the solar constant (the amount of solar energy 
incident at the top of the atmosphere), the factor of 1/4 accounts for the ratio of 
the area the disk on which solar energy is absorbed to the Earth's surface, a = 0.3 
is the albedo (the fraction of solar radiation reflected) . We can solve this equation 
for the "effective temperature" of Earth, T e = 255 K. The surface is much warmer 
than this, as it is heated by the greenhouse effect. 

The general principle of the greenhouse effect is this: the atmosphere de- 
creases in temperature with height, so we can consider the atmosphere to be a 
layer which is colder than surface (T atm < T sur f). This will absorb a fraction of the 
upward travelling radiation with from the surface (the fraction is the absorptivity) . 
Recall that emissivity equals absorptivity. The atmospheric layer emits both up- 
wards and downwards, dependant on its temperature (from the Stefan-Boltzmann 
law) and its emissivity, with = = eaT^ tm . Given that T atm < T sur f, less en- 
ergy is then emitted to space than would have been with no absorbing layer. Energy 
radiated downwards heats the surface. The strength of the greenhouse effect not 
only depends on the amount of greenhouse gas, but the temperature of the layer 
of the atmosphere from which emission occurs. The colder the emitting layer, the 
better the greenhouse effect. 

For the discussion of the runaway greenhouse, some more technical details are 
necessary. A key quantity is the optical depth, t, of a layer or the atmosphere. 
This depends on the number of molecules in the layer on the path the radiation 
takes, N, and some absorption coefficient for the relevant absorbing gas, k, thus 
t = kN. The amount of radiation which is transmitted through a layer depends 
directly on this optical depth: F\ = Foe~ T (Beer's law), where Fo and F\ are 
the incident and transmitted fluxes respectively. Hence the amount of absorption 
depends on the logarithm of the amount of absorbing gas. Further, most of the 
absorption will happen around within an order of magnitude of optical depth of 
one (unity): e" ' 1 = 0.90, c" 1 = 0.37, e~ 10 = 0.000045. 

It is conventional to measure the optical depth of from the top of the atmosphere, 
with r = at the top of the atmosphere, so most radiation is then emitted to space 
from around an optical depth of unity. The strength of the greenhouse effect thus 
depends on the temperature structure of the atmosphere around this level. 

We need to distinguish between condensible and non-condensible green- 
house gasses. Most of the greenhouse gases in Earth's atmosphere (carbon dioxide 
is the dominant non-condensible greenhouse gas) require a much colder temperature 
than is ever experienced in order to condense: we refer to these as non-condensible. 
Water vapour, by contrast, does condense at ambient pressures so is referred to as 
a condensible greenhouse gas. There is a large reservoir of the condensed phase in 
the oceans. The pressure of vapour which is in equilibrium (no net evaporation or 
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condensation) is the saturation vapour pressure. This depends exponentially 
on the ambient temperature. Warming Earth's surface will cause there to be much 
more water vapour in the atmosphere. (It is only because of the temperature of 
Earth's atmosphere that particular gases are condensible or non-condcnsiblc. On 
Mars, which is colder than Earth, carbon dioxide is condensible. On Titan, a moon 
of Saturn, which is colder still, water behaves more like rock docs on Earth and 
methane is a condensible gas.). 

The structure of the atmosphere is important. We genrally consider an atmo- 
sphere with two vertical regions. The lower region, the troposphere, is charac- 
terised by large scale vertical (convective) motions. Given the dominance of convec- 
tion, the vertical temperature profile is set by the behaviour of rising air, following 
a moist adiabatic lapse rate (air cools as it rises, but part of this cooling is 
offset by latent heat release as water condenses from the rising air). In the upper 
region, the stratosphere, convection is negligible and the temperature structure is 
determined solely by the balance of absorbed and emitted radiation; it is said to be 
in radiative equilibrium. The boundary between these layers is the tropopause. 
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